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Using a recently proposed directed random walk test, we systematically investigate the popular 
random number generator RANLUX developed by Liischer and implemented by James. We confirm 
the good quality of this generator with the recommended luxury level. At a smaller luxury level 
(for instance equal to 1) resonances are observed in the random walk test. We also find that the 
lagged Fibonacci and Subtract- with-Carry recipes exhibit similar failures in the random walk test. 

... A revised analysis of the corresponding dynamical systems leads to the observation of resonances in 

00 ' the eigenvalues of Jacobi matrix. 
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I. INTRODUCTION 



Large-scale Monte Carlo simulations need good quality random numbers [y. The half-century long history of 
computer simulations shows successes and failures of many algorithms for pseudo-random number (PRN) generation 

J> I Usually, whenever a new test of randomness is proposed some algorithm for generating PRN's becomes obsolete. 
' A somewhat different point of view is advocated here: using a recently proposed test we argue that all algorithms 
using feedbacks (or lags) have more or less the same defects and belong to the same " universality class of badness" . 
This class of generators, however, should not be rejected and rather improved versions of them can still be considered 
reasonably safe (see the final section). 
QQ , One of the most widely used PRN generators of this class is the shift register (SR) generator ||l2| . Correlations 
^\ • in SR generators were pointed out by Compagner jj] but the warning was ignored by computational physicists who 
^ , usually stick to very practical recipes ||^. 

. It has been observed later by Ferrenberg, Landau and Wong 1^ that the statistical simulations of 2D Ising model 
I ■ by the Wolff single cluster algorithm are very sensitive to the defects of the Kirkpatrick-StoU generator |^ (the SR 
Ph! generator with {p, q) = (250, 103)). 

^ ' Motivated by this difficulty, Liischer Q and James Q have developed a random number generator called RANLUX 
which is based on the subtract-with-carry (SWC) recipe advocated by Marsaglia and Zaman ||^. 

We review here some properties of the generator RANLUX using the random walk test recently developed in |l^] 
by Blote, Heringa and one of the authors. 
■ First, we find analytically that the deviations computed using the random walk test are the same for SWC and 
d ^ the lagged Fibonacci recipes. This finding is in contrast to the common belief that the SWC generator gives better 
quality random numbers than recipes which use the exclusive OR operation (like the SR generators) or addition (like 
the lagged Fibonacci generator). Note, that the RANLUX generator with parameter of luxury equal to is equivalent 
to the SWC generator. 

Second, we have found numerically, that correlations among random numbers on distances associated with the lag 
values are still visible for luxury level parameter equal to 1. However at luxury level equal to 2 the correlations lie on 
the boundary of observed deviations (and this is in agreement with the Liischer and James results). 

Therefore, our third result, supporting the author's expectations, is that the RANLUX generator with luxury levels 
higher than 2 is still safe in runs using less than 10^^ random numbers. 

Next, we also revise the dynamical system description given by Liischer in his original paper. Namely, we argue 
that the correlations responsible for the deviations found in the simulation of the two-dimensional Ising model and 
of the random walks are not due to the correlation of trajectories in the corresponding dynamical systems. These 
deviations are rather connected with the correlations between random numbers, on distances equal to the lag values, 
as we show using the random walk test. 

We compute the full spectrum of eigenvalues for the RANLUX generator and for the lagged Fibonacci generator 
with different lags. Our fifth result is the observation of resonances in the eigenvalues, reflecting the complexity of 
the phase space in the area preserving (Hamiltonian) dynamical systems. 
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The paper is organized is follows. In section |l| we recall the definitions of random number generators used here 
- Subtract-with-Carry (SWC), lagged Fibonacci (LF) and Shift register (SR), and review the random walk test. In 
section III we apply the random walk test to the RANLUX generator and surprisingly find that the SWC generator 



obeys the same analytical solutions as the LF generator. Then we analyse how the correlations in random walks 
depend on the luxury levels of the RANLUX generator. In section |^ we reanalyse the dynamical system approach to 
RANLUX generator and find resonances in the eigenvalues of the Jacobi matrix associated with the complex structure 
of the manifolds. Finally, in section^ we discuss the attempts to improve the lagged random number generators. 



II. DEFINITIONS AND NOTATIONS 



We recall here the definitions of several random number generators which have the common feature of using lag 
algorithms. We also recall the definition of the random walk test. 



A. Marsaglia-Zaman recipe and RANLUX generator 



The PRN generation algorithm of Marsaglia-Zaman is defined by a recursion relation involving three fixed 
positive integers 6, r, s where b is called the base and r > s are called the lags. This algorithm is usually known as 
SWC generator. Given the first r PRN's xo,xi, ...Xr-i and the "carry bit" c,_i, the n-th PRN (n > r) is given by 

Xji — Xyi^j- Cyi— i) mod b 

Cn — ■^71 — s -^n—r ^0 /-< \ 

and 

c„ 1 otherwise 



The maximum possible period of this generator is M — b^ — b^ + 1 and it is attained when M is prime and 6 is a 
primitive root modulo M. 

In order to improve the properties of this algorithm which fails to pass some correlation tests, Liischer proposed 
to discard some of the PRN's produced by the MZ recursion and to use only the remaining ones. James has then 
defined ||^ four levels of rejection (called " luxury levels" ) , characterized by an integer p > 24 in which the generator 
produces 24 PRN's, then discards the successive p — 24 and so on. Clearly the value p = 24 reproduces the original 
MZ recipe where all PRN are kept and it is called luxury level zero. The values p = 48, 97, 223, 389 define the luxury 
levels 1, 2, 3, 4 respectively. It has been suggested that level 3 has a good chance of being optimal. 



B. Lagged Fibonacci generator 



The Fibonacci lagged generator (LF) is defined by the recursion relation 

•^n — i'^n — s ^ '^n — r^ mod 2 



(2) 



where again s and r > s are the lags and 2™ is the base. Brent gave conditions |11|, under which the above relation 
generates a sequence of PRN's having the maximum possible period T — 2™~^(2'' — 1). (T is simply the product of 
the reduced computer word length and of the number of possible initial seeds. The latter is the number of possible r 
bit seeds excluding zero.) 



C. Shift Register generator 

The generalized feedback shift register (SR) method is defined by the recursion rule 

•^n — '^n — s ® •^n — r (3) 

where ® denotes the exclusive-or operation (bitwise sum modulo 2), s and r > s are the lags. This sequence has 
maximum period T = 2*" — 1 if its characteristic polynomial is irreducible (for details, see Chapter II of Golomb book 
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D. Random Walk Test 



Let us consider the one-dimensional directed random walk model |10|: a walker starts at some site of an one- 
dimensional lattice and, at discrete times i, either he takes a step in a fixed direction with a probability fii or he stops 
with a probability 1 — fii. In the latter case a new walk begins. The probability of a walk with length n is then 



In the case in which all probabilities are equal = ^, the the probability for a walk of length n is 

P{n) = (1 - /.) (4) 

and the mean walk length is given by 

, . 1 



(5) 



Comparing /i with tanh{J/kBT) it immediately follows that the mean walk size diverges as /i — > 1 in the same way 
as the mean cluster size diverges as the temperature T ^ in the ID Ising model pM. 



III. RANDOM WALK IN A CORRELATED ENVIRONMENT 



The production rules for the the SR and LF generators with the pair {p, q) of lags give rise to correlations on 
distances Ic — ip + jq, where i and j are integers P,p^. We can treat these correlations as some kind of impurities, 
placed on the bonds of ID lattice and thus influencing the walk probabilities fii. 

In some cases, it is possible to get exact analytical results for the deviations of walk probabilities from their values 
in the pure case. 

In addition, correlations in the PRN sequences can be nicely detected by simulation of this model on the computer 
and by comparing the observed frequency P*{n) of walks of length n with their expected (pure) probability P{n). In 
what follows we will plot the measured deviations 

5Pin)^^-l (6) 
P(n) 

versus n, for various PRN generators and briefly comment on why they occur. 

For all PRN generators we will consider here, biases are brought out in the directed random walk model by the 
obvious fact that the end of a walk and the beginning of a new one are strictly correlated: as a result probability 
deviations will occur at walk lengths related to the lag values in the production rule. To be definite, suppose for 
example that in our PRN sequence the n-th term x„ is generated from Xn-r and a;„_s with r > s and let us assume 
for simplicity that the stepping probability in our model has the value /i = 1/2. In other words the ideal walker would 
proceed after tossing a good coin, while the real walker will have to resort to our deterministic PRN generator with 
lags r and s. Let us now suppose that the real walker stops after k steps, so that Xk > /i. Since the production rule 
gives Xk in terms of Xk-r and x^-s, both of which have to be less than /x, it follows that the probability of the k-th 
step of the real walker is biased by his previous history. 

The first study of this new correlation test in the case of the shift register and lagged Fibonacci generators has been 
performed in pO[ |. 



A. Random walk and luxury levels 



Here we shall test the RANLUX generator with different levels of luxury in the simulation of the directed random 
walk model. 

To begin let us choose luxury level and step probability /i = 31/32. In Fig. 1 we have plotted versus n the 
deviation SP{n) defined by of the frequency of the walks of length n from their expected probability. We notice 
strong signals for ni = r = 24 and n2 = 2r = 48. A sizable deviation is also clearly visible for all walk lengths 
exceeding the largest lag r. 
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We should stress that our choice of the value of the step probability has no special meaning and that our results 
depend smoothly on fi. 

One should notice the qualitative similarity of the deviations in Fig. |l| and Fig. ^, which refers to the case of a 
lagged Fibonacci generator. Moreover, it could be shown analytically that the deviations aX n — r and at n — i, 
satisfying 

( r + s a s < ^ 

r <i< < -f ^ ^ (7) 

1 2r — s it s > ^ ^ ' 

coincide for the lagged Fibonacci and the SWC generators. We recall that the RANLUX generator with luxury level 
is equivalent to the SWC generator. It also could be shown, that the mean value of the deviation for the SWC 



generators at n = r is equal to that for the lagged Fibonacci generator obtained in 1 10 

r 



SPi-) - ^ (8) 



and that the mean value of the deviation at n = r + 1 is 

5F(r + l) = i^^i-^-l. (9) 

In order to check that point, in Table | we have presented a comparison of the observed values of the deviations 
obtained using either the RANLUX generator with luxury level or the lagged Fibonacci generator for two values of 
the step probability with the expected values given by (|) and (|). All data agree within the statistical errors. 

In contrast to common belief this result shows that there are no differences in the accuracy of the data obtained 
using the LF and the SWC generators, at least as far as the cluster formation processes or the random walks are 
simulated. 

The next luxury level of Liischer and James RANLUX is 1 . In this case only the first r — 2A random numbers are 
used out of each generated set of p = 48. The observed deviations are plotted in Fig. ||. Clearly, a deviation at 
n = r — 24 still remains, but has opposite sign and a much smaller size in comparison with that observed for the 
luxury level 0. 

Also at the values n — 28, 34, 38, 48, all of which are linear combinations of the lags with small integer coefficients, 
smaller deviations remain clearly visible. 

We conclude that, at the luxury levels and 1, the deviations have the same qualitative nature, the only difference 
being that the absolute value of the deviations is depressed for the higher levels of luxury. 

This remark is confirmed by Fig. ^ where we have plotted the observed deviations for the case of luxury level 2, and 
the visible deviations are not larger than the statistical errors. Notice that in this case the average length of the walk 
is 32 and for the total number of 10^^ walks approximately 10^^ random numbers were used. We could extrapolate 
our result in like in Ref. and we expect that deviations will be visible even for the luxury level 2 if « 10^^ random 
numbers are used. 



B. Lagged Fibonacci with luxuries 

In Fig. I we present the results of a simulation using the plain lagged Fibonacci generator with lags r = 24 and 
s = 10. The plot shows features completely similar to those of Fig. |l| obtained with RANLUX at luxury level 0. 

We can now extend to the lagged Fibonacci generator (LF) the procedure of decorrelation which discards 24 PRN's 
out of each generated set of 48, and call it following Liischer Q and James |^ luxury level 1. In Fig. || we have 
reported the results of a simulation using the LF generator with lags r = 24, s = 10 at luxury level 1. The results 
again show no qualitative difference from those reported in Fig. ^ obtained when RANLUX at the same luxury level 
was used: a strong deviation is observed when the walk length equals the value of the largest lag r and smaller but 
visible deviations occur for any n > r. 

In order to understand how the decorrelation procedure works we also have studied intermediate luxury levels 
obtained by varying p between 24 and 48 in steps of one. In Fig. ^ we have plotted 5P(24) versus p — r. The value 
of SP{2A) decreases rapidly but peaks are still evident at p — r = 10, 14, 24 indicating the persistence of an interplay 
between p and the lags r and s. 
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IV. DYNAMICAL SYSTEM 



Liischer justified his proposal of skipping some p — r numbers using the language of dynamical systems 
Let us then briefly recall here a few well known results from the dynamical system theory. 
The production rules (|l|),(|^) and (||) can be rewritten in the matrix form 



L v„, (10) 



where the rule of matrix "multiplication" includes the "®" operation for the shift register generator, the 
usual "+" operation for the lagged Fibonacci and the "— " operation for the SWC generators, and v„ = 
{Xn-r+i, ■■■,Xn-i,Xn) is an T-component vector. 

The simplest example is the Fibonacci sequence 

Xn ^ Xn-l + (11) 



which can be written in matrix form 



Xn ] = ( ^ ^] ( ^"-1 
Xn-1 1 V 1 / V 



(12) 



This shows that the dynamical system, corresponding to the Fibonacci random number generator with lags r — 2 and 
s = 1 is the famous Arnold's "cat map", acting on the unit torus, obtained by rewriting (Kw in the form 



v„+2 = M v„. (13) 
The matrix = L*" is called the monodromy matrix, 

The eigenvectors (7"^,1) and (7^ — 1,1) of M define respectively the stable and the unstable directions of the 
origin, which is an unstable fixed point. The corresponding eigenvalues are 



Ai \ / 72 \ / 2.618034.. 
A2 j ^ I 7-2 ) ~ I .381966... 



(15) 



and 7 = (1 + \/5)/2 is the golden mean. 

Fig. 1^ demonstrates the idea of hyperbolicity of the "cat map". The points on the part / of the unstable manifold 
are mapped under a single application of the transformation ([l^) onto the four parts (in order, denoted by I, II, III 
and IV) of the same unstable manifold. Clearly, Fig. ^ shows a trajectory winding on the torus with the golden angle 
tan (j) — {VE — l)/2. The distance between neighbor points along the unstable manifold is enlarged by a factor Ai. 
In order to imagine the complexity of the trajectories, one needs to add the set of transversal lines representing the 
stable manifold. Along these lines the distances are contracted by a factor A2 . Thus, according to the Poincare theory, 
there is an infinite set of periodic points with zero measure in this system (we address interested readers to the figures 
on page 170 of Ref. 0.) 

It is known [l^ that the "cat map" is an Anosov system, which implies the following properties of randomness 
|p7| : global instability (the Lyapunov exponent, which is the logarithm of Ai from Eq. |l^, is positive), positivity 
of Kolmogorov-Sinai entropy (the trajectories are locally divergent), mixing (existence of equilibrium state for the 
distribution function) and, finally, ergodicity (the equivalence of averaging over time and over space). Nevertheless, 
the "cat map" doesn't have the strongest property of randomness, namely, it is not a Bernoulli system, in other words 
(M. Tabor in Ref. [jl^, page 174) it is not a system "whose motion is as random as a fair coin toss". 

In fact, the lagged Fibonacci generators (and the Marsaglia-Zaman SWC generators as well) are an r-dimensional 
generalization of a "cat map". A regular way to extend randomness properties to higher dimensions is not known. 
Nevertheless, in general, one should not expect from the generalized system stronger properties of randomness than in 
the case of a low dimensional system. In fact, practically it is sufficient to analyze the eigenvalues of the corresponding 
monodromy matrix to be sure that a system is locally hyperbolic in any point of phase space (to have, at least, the 
properties of the weakly Anosov systems jl^). 

The suggestion of Liischer to improve the quality of the random numbers generated by rule ( |l0[ ) is based on the 
observation, that since the absolute value of the largest eigenvalue of the matrix L in (HO) is greater than unity, we can 



5 



construct a monodromy matrix A4 = li^ with a relatively large modulus of the largest eigenvalue, whose logarithm is 
the Lyapunov exponent of the corresponding dynamical system on the r-dimensional torus. 

Indeed, his Fig. 1 from Ref. ||^ shows clearly the process of decorrelation (loos of memory). If we apply the 
monodromy matrix M to the shortest representable vector which has elements of order of 2^™, the length of our 
vector will become of order of unity after a number of iterations n such that 

2"™ A" w 1. (16) 

Knowing Liischer result for n = 16 and considering the case of 24-bit representation for the random numbers (to = 24), 
as in his paper, from ( |l6| ) we could estimate A as A « 2^^/^^ « 2.828... in a good agreement with his direct calculation 
A = 2.746.... 

This means, by Liischer, that the decorrelation time depends on the number of bits in the random number repre- 
sentation 

n^TOj^^ (17) 
log A 

and leads to the prediction of an infinite decorrelation time in the limit to — s- oo. We could use the same arguments 
for the case of a "cat map" and would arrive to a contradiction with Arnold's result, that the "cat map" is a C-system 
with the strong ergodic properties cited above. 

This means, that the values of the Lyapunov exponent are not directly connected with the process of decorrelation 
in random number sequences. We can support our idea considering LCG generator which are known to have a lattice 
structure and introducing luxury levels in this case. It is clear, that the lattice structure will not disappear - increasing 
the luxury levels would enlarge eigenvalues, but the lattice structure would still persist. 

Really, as we have demonstrated using the random walk test, the situation is more complicated and we need to 
consider also the correlations in the production rules on distances equal to the lags and to their linear combinations. 
This effect is clearly visible in Fig. ^, where one can see strong resonances at lag values r = 24 and s = 10, whereas 
the eigenvalues are increasing smoothly as Af with p. 

It should be noted, that the deviations at dP(r) do not vanish smoothly on the way from one luxury level to another, 
but they oscillate. We hope, nevertheless, that the amplitude of these oscillations vanishes. 



A. Eigenvalue spectrum 

Each state of the random number generators of interest is described by a set of r random numbers , (i = 1, 2, r). 
After r applications of any of the rules (|l|),(||) and (H), we will have completely refreshed the set of r random number 
u^, (i = 1, 2, r). This defines the mapping of unit r-dimensional cube onto itself 

We could construct numerically the Jacobi matrix p6|jT9[ , as the matrix of the variations of the final set Ui = u^+Sui 
due to the variations of the initial set Vi = + Svi: 



J., = ^. (18) 



which reflects the behavior of neighbor trajectories for all components of the set Ui. 

The eigenvalues of RANLUX with luxury level are plotted on Fig ^ together with the eigenvalues for corresponding 
lagged Fibonacci generator with lags (24, 10). It is clear that eigenvalues will collapse on a single curve changing the 
sign of i?e(A). These eigenvalues could be determined directly from the characteristic polynomial for the productions 
rules, which reads as 

A24-Ai4±l = (19) 

with plus sign for the Marsaglia-Zaman and minus sign for the Fibonacci generator. In Fig. ^ we have plotted A^** 
(we remind, that the dynamical system is obtained as the mapping of the r-dimensional unit cube onto itself, after p 
iterations of the production rule. The luxury level corresponds to p = 24.) 

Fig. ^ demonstrates how the eigenvalues change going from the luxury level to the luxury level 1 in the case of 
the RANLUX generator. 

Below we will produce data for the lagged Fibonacci generator with p > r (we generate p random numbers and use 
only the last r of them) taking into account the similarity in many aspects of the lagged Fibonacci and the RANLUX 
generators described in this paper. 
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First, we show how the eigenvalues depend on the value of the lag r. Fig. |lO| shows, in particular, that the eigenvalues 
for (24,10) and (250,103) lie almost on the same curve. We will plot below the eigenvalues for the last case in order 
to have more detailed pictures. 

In Figs. pl - [r^ we show how the eigenvalues change when the parameter p is varied. One clearly sees that the 
resonances in eigenvalues reflect the interplay of the lags (r, s) and the parameter p. These pictures reflect in some 
sense the complexity of the phase space of our dynamical systems. 



V. CONCLUSIONS 



In this paper we have shown that the sequences of random numbers generated using both the lagged Fibonacci and 
the subtract-with-carry generators have equivalent correlation properties, at least as far as the random walk and the 
cluster algorithms are concerned. 

We have argued that correlations in random numbers on the lag distances decrease with increase of the luxury levels, 
although not monotonically. There are also resonances when the number p — r of discarded random numbers coincides 
with linear combinations of lags. We also have argued that correlations on these distances are more important than 
correlations between "trajectories" of random numbers in phase space. 

Comparisons of lagged Fibonacci and shift register generators have been presented in pO[ | and our analysis demon- 
strates that they behave in a qualitatively similar way. 

One can conclude, that all recipes using lags will produce sequences of random numbers with the same qualitative 
features. It is important that knowing these bad properties we can avoid possible deviations either using the idea of 



luxury levels [Qj^li or decimated sequences or combinations of two or more random number generators 14 

or generators with a larger number of lags P,]l5| . 

Therefore rephrasing Orwell: "All algorithms for generating random numbers are equally bad, but some of them 
are more equal than others" pO|). 
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TABLE L Comparison of probability deviations, computed using directed random walk test with lagged Fibonacci (LF) and 
subtract-with-carry (SWC) random number generators, and predicted by formulas ^ and |^. 



1 A* = 31/32 


At = 15/16 




5P{r) 


5P{r + 1) 


5P{r) 


5P{r+l) 


Expected 


-0.48387... 


0.03146... 


-0.4666... 


0.06319... 


LF 


-0.48304(55) 


0.03221(81) 


-0.46636(47) 


0.06290(85) 


SWC 


-0.48299(64) 


0.03122(81) 


-0.46696(63) 


0.06243(113) 
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FIG. 1. Deviation 5P of the probability of a walk length n from the value for uncorrelated random numbers versus walk 
length. We have used the RANLUX random number generator with lag values r = 24 and s = 10 and luxury level 0. The step 
probability is /u = 31/32 = 0.968750. The result of averaging over 10® walks is shown. 
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FIG. 2. Deviation SP of the probability of a walk length n from the value for uncorrelated random numbers versus walk 
length. We have used the RANLUX random number generator with lag values r = 24 and s = 10 and luxury level 1. The step 
probability is // = 31/32 = 0.968750. The result of averaging over 10^^ walks is shown. 



9 



0.0005 
0.0004 
0.0003 - 
0.0002 
0.0001 

': 

-0.0001 
■0.0002 

-0.0003 
-0.0004 - 
•0.0005 



5 10 15 20 25 



35 40 45 50 55 60 



FIG. 3. Deviation SP of the probability of a walk length n from the value for uncorrelated random numbers versus walk 
length. We have used the RANLUX random number generator with lag values r = 24 and s = 10 and luxury level 2. The step 
probability is = 31/32 = 0.968750. The result of averaging over lO'^^ walks is shown. 
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FIG. 4. Deviation dP of the probability of a walk length n from the value for uncorrelated random numbers versus walk length. 
We have used the lagged Fibonacci generator with lag values r = 24 and s = 10. The step probability is /j, = 31/32 = 0.968750. 
The result of averaging over 10* walks is shown. 
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FIG. 5. Deviation 5P of the probability of a walk length n from the value for uncorrelated random numbers versus walk 
length. We have used the lagged Fibonacci generator with lag values r = 24 and s = 10 and luxury level 1. The step probability 
is /U = 31/32 = 0.968750. The result of averaging over 10^° walks is shown. 
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FIG. 6. Lagged Fibonacci generator with lags r = 24 and s = 10 and ^ = 15/16. The figure shows the deviation 6P{r) as a 
function of p — r on the way from luxury level (p — r = 0) to luxury level 1 {p — r = 24). 
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FIG. 7. Fibonacci sequence on torus: Arnold's cat map. 
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FIG. 8. Eigenvalues of the Jacobi matrixes corresponding to the RANLUX generator with luxury level (marked with pluses 
+). The same for lagged Fibonacci generator with lags (24,10) (marked with diamonds O). 
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FIG. 9. Eigenvalues of the Jacobi matrixes corresponding to the RANLUX generator with luxury level (pluses +) and to 
the RANLUX generator with luxury level 1 (diamonds O). 
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FIG. 10. The complex plane eigenvalues of the Jacobi matrices corresponding to the lagged Fibonacci generator with lags 
(24,10), (36,11), (89,38), (127,64) and (250,103) indicated by (□,-!-, x,*, □) respectively. 
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FIG. 11. Complex plane eigenvalues for the lagged Fibonacci generator with lags (250,103). The number of discarded PRN's 
is 3, 5 and 22 going from top to bottom. 
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FIG. 12. Complex plane of eigenvalues for the lagged Fibonacci generator with lags (250,103). The number of discarded 
PRN's is 49, 50 and 54 going from top to bottom. 
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FIG. 13. Complex plane of eigenvalues for the lagged Fibonacci generator with lags (250,103). Number of discarded PRN's 
is 71, 92 and 125 going from top to bottom. 
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FIG. 14. Complex plane of eigenvalues for the lagged Fibonacci generator with lags (250,103). The number of discarded 
PRN's is 152, 184 and 206 going from top to bottom. 
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